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Abstract 

The nonleptonic heavy meson decays B ^ *^7r(p), and D —^ 

Ki*) TV are studied based on the three-scale perturbative QCD factorization 
theorem developed recently. In this formalism the Bauer-Stech-Wirbel pa¬ 
rameters Oi and 02 are treated as the Wilson coefficients, whose evolution 
from the W boson mass down to the characteristic scale of the decay processes 
is determined by effective held theory. The evolution from the characteristic 
scale to a lower hadronic scale is formulated by the Sudakov resummation. 
The scale-setting ambiguity, which exists in the conventional approach to 
nonleptonic heavy meson decays, is moderated. Nonfactorizable and non¬ 
spectator contributions are taken into account as part of the hard decay 
subamplitudes. Our formalism is applicable to both bottom and charm de¬ 
cays, and predictions, including those for the ratios R and Rl associated 
with the B —>■ J/ip decays, are consistent with experimental data. 
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I. INTRODUCTION 


The analysis of exclusive nonleptonic heavy meson decays has been a chal¬ 
lenging subject because of the involved complicated QCD dynamics. These 
decays occur through the Hamiltonian 

H = , ( 1 ) 

with Gp the Fermi coupling constant, U’s the Cabibbo-Kabayashi-Maskawa 
(CKM) matrix elements, g’s the relevant quarks and (qq) = ^7/^(1 — J 5 )q the 
V — A current. Hard gluon corrections cause an operator mixing, and their 
renormalization-group (RG) summation leads to the effective Hamiltonian 

H.S = ^l/yVS|ci(f,)Oi + C2(f,)02| , (2) 

where the four-fermion operators Oi ,2 are written as Oi = {qiqk){qjqi) and 
O 2 = {Qj(lk){qiqi)- The Wilson coefficients ci, 2 , organizing the large loga¬ 
rithms from the hard gluon corrections to all orders, describe the evolution 
from the W boson mass Mw to a lower scale /i with the initial conditions 
ci{Mw) = 1 and C 2 {Mw) = 0. 

The simplest and most widely adopted approach to exclusive nonleptonic 
heavy meson decays is the Bauer-Stech-Wirbel (BSW) model |l[ based on 
the factorization hypothesis, in which the decay rates are expressed in terms 
of various hadronic transition form factors. Employing the Fierz transfor¬ 
mation, the coefficient of the form factors corresponding to the external W 
boson emission is oi = ci -1- C 2 /N, and that corresponding to the internal W 
boson emission is a 2 = C 2 + Ci/N, N being the number of colors. The form 
factors may be related to each other by heavy quark symmetry, and be mod¬ 
elled by different ansatz. The nonfactorizable contributions which can not be 
expressed in terms of hadronic transition form factors, and the nonspectator 
contributions from the W boson exchange (or annihilation) are neglected. In 
this way the BSW method avoids the complicated QCD dynamics. 

Though the BSW model is simple and gives predictions in fair agreement 
with experimental data, it encounters several difficulties. It has been known 
that the large N limit of 01 ^ 2 , the choice ai = Ci(Mc) ~ 1.26 and 02 = 
C 2 {Mc) ~ —0.52, with Me the c quark mass, explains the data of charm 
decays |]^. However, the same large N limit of oi = ci(Mfe) Ri 1.12 and 
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02 = C 2 {Mb) ~ —0.26, Mb being the b quark mass, does not apply to the 
bottom case. Even after including the 01 ^ 2 /N term so that oi = 1.03 and 
02 = 0.11, the BSW predictions are still insufficient to match the data. To 
overcome this difficulty, parameters y, denoting the corrections from the 
nonfactorizable hnal-state interactions, have been introduced p|. They lead 
to the effective coefficients 



(3) 


X should be negative for charm decays, canceling the color-suppressed term 
1/N, and be positive for bottom decays in order to enhance the predictions. 
Unfortunately, the mechanism responsible for this sign change has not been 
understood completely. Furthermore, in such a framework theoretical predic¬ 
tions depend sensitively on the choice of the scale fi for the Wilson coefficients: 
Setting fi to 2Mb or Mb/2 gives rise to a more than 20% difference. 

Equivalently, one may regard 01^2 as free parameters, and determine them 
by data htting. However, the behavior of the transition form factors involved 
in nonleptonic heavy meson decays requires an ansatz |Q as stated above, 
such that the extraction of 01,2 from experimental data becomes model- 
dependent. On the other hand, it was found that the ratio 02/01 from an 
individual £t to the CLEO data of H ^ Zl^*^7r(p) P varies significantly 
d]. It was also shown that an allowed domain ( 01 , 02 ) exists for the three 
classes of decays —> Zi)h)+^ ^0 and B~ —> Zi)h)o^ when the 

experimental errors are expanded to a large extent 0. 

Moreover, it has been very difficult to explain the two ratios associated 
with the B —>• J/'ipK^*^ decays Q], 

B{B ^ J/^PK*) 

B{B ^ J/'i/K) ’ ^ ^ Jl-i/K*) ’ ^ ’ 

simultaneously in the BSW framework, where B{B —> J/i/K) is the branch¬ 
ing ratio of the decay B —> J/'ipK. It was argued that the inclusion of nonfac¬ 
torizable contributions is essential for the resolution of this controversy [0]. 
Such contributions have been analyzed in 0, || based on the Brodsky-Lepage 
approach to exclusive processes IQ, in which the full Hamiltonian in Eq. (0) 
was employed. It was found that the nonfactorizable internal lU-emission 
amplitudes are of the same order as the factorizable ones. However, this 
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naive perturbative QCD (PQCD) approach can not account for the destruc¬ 
tive interference between the external and internal hh-emission contributions 
in charm decays. This is obvious from the fact that the coefficient associated 
with the internal W emissions is 02 = 1/A^ in both bottom and charm decays, 
and thus does not change sign. 

Recently, a modified PQCD formalism has been proposed following the 
series of works 0,0, [1,111, im, where the PQCD formalism constructed 
from the full Hamiltonian H was shown to be applicable to the B ^ de¬ 
cays 10] in the fast recoil region of final-state hadrons [0. It was recognized 
that nonleptonic heavy meson decays involve three scales: the W boson mass 
Mw, the typical scale t of the decay processes, and the hadronic scale of or¬ 
der Aqcd- Accordingly, the decay rates are factorized into three convolution 
factors: the “harder” hP-emission function, the hard b quark decay subampli- 
tude, and the nonperturbative meson wave function, which are characterized 
by Mw, t and Aqcd, respectively. Radiative corrections then produce two 
types of large logarithms ln(Mvc/t) and ln(t/AQCD)- In this three-scale fac¬ 
torization theorem \a{Mw/t) are summed to give the evolution from Mw 
down to t described by the Wilson coefficients ai, 2 (t), and ln(t/AQCD) are 
summed into a Sudakov factor 0, which describes the evolution from t 
to the lower hadronic scale. The former has been derived in effective field 
theory, and the latter has been implemented by the resummation technique 

H- 

This modified PQCD formalism is p-independent, ie. RG-invariant, and 
thus the scale-setting ambiguity existing in conventional effective field the¬ 
ory is moderated 0|. As the variable t runs to below Mi, and M^, the con¬ 


structive and destructive interferences involved in bottom and charm decays, 
respectively, appear naturally. Furthermore, not only the factorizable, but 
the nonfactorizable and nonspectator contributions are taken into account 
and evaluated in a systematic way. With the inclusion of the nonfactorizable 
contributions, we find that 01^2 restore their original role of the Wilson coef¬ 
ficients, instead of being treated as the BSW free parameters. The branching 
ratios of various decay modes B ^ *^7r(p) and D —^ TT, and the ratios 
R and Rl associated with the B J/ipK^*'^ decays can all be well explained 
by our formalism. 

In Sec. II we derive the three-scale PQCD factorization theorem, con¬ 
centrating on the separation of the contributions characterized by different 
scales. The incorporation of the Sudakov resummation is briefly reviewed. 
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In Sec. Ill the decays B ^ 7r(p) are investigated to demonstrate the im¬ 

portance of the nonfactorizable contribntions. We then apply the formalism 
to the decays D in Sec. IV, and show that the internal hh-emission 

amplitnde can become snfliciently negative in charm decays. In Sec. V we 
compnte the decay rates of i? —>■ J/'ipK^*'^ and find that the predictions for 
R and Rl match the data simultaneously. Section VI is the conclusion, 
where possible further improvements and applications of our approach are 
proposed. 


II. THREE-SCALE FACTORIZATION THEOREMS 


In this section we construct the modified PQCD formalism, that embod¬ 
ies both effective field theory and factorization theorems. The motivation 
comes from the fact that the Wilson coefficients ci ,2 of the effective Hamilto¬ 
nian in Eq. (H) are explicitly p-dependent. Since physical quantities such as 
the decay rates, which are expressed as the products of ci ,2 with the matrix 
elements of the four-fermion operators Oi^ 2 , do not depend on /i, the latter 
should contain a p dependence to cancel that of the former. However, such a 
cancellation has never been implemented in any previous analysis of nonlep- 
tonic heavy meson decays. As stated in the Introduction, the BSW method 
employs the factorization hypothesis [Q, under which the matrix elements of 
Oi 2 are factorized into two hadronic matrix elements of the (axial) vector 
currents {qq). Since the current is conserved, the hadronic matrix elements 
have no anomalous scale dependence, and thus the /i dependence of the Wil¬ 
son coefficients remains. To remedy this problem, /i should be chosen in 
such a way that the factorization hypothesis gives dominant contributions. 
However, the hadronic matrix elements involve both a short-distance scale 
associated with the heavy quark and a long-distance scale with the mesons. 
Naively setting /i to the heavy quark mass will lose large logarithms contain¬ 
ing the small scale. It is then quite natural that theoretical predictions are 
sensitive to the value of p [ 1^5 HI • 

We shall show that the cancellation of the /i dependence is explicit in our 
formalism. We begin with the idea of the conventional PQCD factorization 
theorem for the B —>■ D transition form factors, which describe the amplitude 
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of a 6 quark decay into a c quark through the current operator {cL’^^hL). Ra¬ 
diative corrections to these form factors are ultraviolet (UV) finite, because 
the current is not renormalized. However, the corrections give rise to infrared 
(IR) divergences at the same time, when the loop gluons are soft or collinear 
to the light partons in the mesons. These IR divergences should be sepa¬ 
rated from the full radiative corrections and grouped into nonperturbative 
soft functions. 

The separation of IR divergences in one of the higher-order diagrams is 
demonstrated by Fig. 1(a), where the bubble represents the lowest-order de¬ 
cay subamplitude of the B meson. This diagram is reexpressed into two 
terms: The first term, with proper eikonal approximation for quark prop¬ 
agators, picks up the IR structure of the full diagram. The second term, 
containing an IR subtraction, is finite. The first term, being IR sensitive, 
is absorbed into a meson wave function 0(6,/i), if the radiative correction is 
two-particle reducible, or into a soft function t/(6,/i), if the radiative correc¬ 
tion is two-particle irreducible as shown in Fig. 1(a). Here b is the conjugate 
variable of the transverse momentum kx carried by a valence quark of the 
meson, and thus can be regarded as the transverse extent of the meson. It 
will become clear later that the scale 1/6 serves as an IR cutoff of the as¬ 
sociated loop integral. The function U corresponds, in some sense, to the 
nonfactorizable final-state interactions in the literature of nonleptonic heavy 
meson decays p[. The second term, being IR finite, is absorbed into the 
hard decay subamplitude as a higher-order correction, where t is its 

typical scale. 

The above factorization procedure is graphically described by Fig. 1(b), 
where the diagrams in the hrst parentheses contribute to H, and that in the 
second parentheses to 0 or U. Below we shall neglect U, because of the pair 
cancellation between the diagram in Fig. 1(a) and the diagram with the right 
end of the gluon attaching the lower quark line. The combination of these 
two diagrams leads to an integrand proportional to a factor 1 — It being 
the transverse loop momentum. It is then apparent that U is unimportant, 
if the main contributions to the form factors came from the small 6 region. It 
will be shown that the Sudakov factor mentioned in the Introduction exhibits 
a strong suppression at large 6, and thus justifies the neglect of U. On the 
other hand, U involves complicated color flows. Hence, we leave its discussion 
to a separate work ii. 

Though the full diagrams are UV hnite, the IR factorization introduces 
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UV divergences into 0 and H, which have opposite signs. This observation 
hints a RG treatment of the factorization formula derived above. Let 7 ,^ be 
the anomalous dimension of 0. Then the anomalous dimension of H must 
be — 7 ^. We have the RG equations 

d d 

/ X —0 = -70 = , 5 

dfi dfi 

whose solutions are given by 


0(6,/i) 


1 / 6 ) exp 
exp 


- f —70(«.(h)) 
Ji/b /i 

70(«s(h)) 


dfi 
'i/b fi 

R dfi 

I fl jJj 


( 6 ) 

m 


Equation (^) describes the evolution of 0 from the IR cutoff 1/6 to an arbi¬ 
trary scale fi, and Eq. m describes the evolution of H from fi to the typical 
scale t. The physics characterized by momenta smaller than 1/6 is absorbed 
into the initial condition 1 / 6 ), which is of nonperturbative origin. After 
the RG treatment, the large logarithms ln(t//x) in H are grouped into the 
exponent, and thus the initial condition H[t, t) can be computed by perturba¬ 
tion theory. Gombining Eqs. (j^) and (|^), the factorization formula becomes 
free of the fi dependence as indicated by 


H{t, fi)(j){b, fi) = H(t, t)(f){b, 1/6) exp 



( 8 ) 


The effective Hamiltonian ifeff in Eq. (|]) can be constructed in a similar 
way. Gonsider the nonleptonic 6 quark decays through a W boson emission 
up to 0{as). We reexpress the full diagram, which does not possess UV 
divergences because of the current conservation and the presence of the W 
boson propagator, into two terms as shown in Fig. 2(a). The hrst term, 
obtained by shrinking the W boson line into a point, corresponds to the 
local four-fermion operators Oi ^2 appearing in Hgs, and is absorbed into the 
hard decay subamplitude H(t,fi). This subamplitude is characterized by 
momenta smaller than the W boson mass that is, by the typical scale t 
of the heavy meson decays, since gluons in H do not ’’see” the W boson. The 
second term, characterized by momenta of order Mw due to the subtraction 
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term, is absorbed into a ’’harder” function Hr{Mw,^^) (not a amplitude). 
Note that the factorization in H is not complete yet, because it still contains 
IR divergences, ie. the contributions characterized by the hadronic scale. 

We then obtain the 0{as) factorization formula shown in Fig. 2(b), where 
the diagrams in the hrst parentheses contribute to H^, and those in the second 
parentheses to H. This formula should be interpreted as a matrix relation 
because of the mixing between the operators Oi and O 2 , or equivalently, 
the four-fermion vertex should be regarded as the linearly combined opera¬ 
tors Oi ± O 2 , which evolve independently. The four-fermion vertex in the 
denominator means that Hr does not carry Dirac and color matrix struc¬ 
tures. Similarly, UV divergences are introduced in the above factorization 
procedure, when the W boson line is shrunk, and thus both H and Hr need 
renormalization. The RG improved factorizaton formula is written as 


Hr{Mw, l^)H{t, jj,) = Hr{Mw, Mw)H{t,t) exp 


dii ■ 

/ —^HriasW) 
Jt fi 


. (9) 


with the anomalous dimension of Hr- It is easy to identify the exponential 
in Eq. (^ as the Wilson coefficient, implying that the scale p. in the Wil¬ 
son coefficient should be set to the hard scale t. The function Hr{Mw, Mw) 
can now be safely taken as its lowest-order expression H^^'^ = 1, since the 
large logarithms ln(Miy//i) have been organized into the exponent. Note 
that the appropriate active flavor number should be substituted into as(/2), 
when /i evolves from Mw down to t. The continuity conditions for the tran¬ 
sition of as{fl) between regions with different active flavor numbers PO are 
understood. 

We are now ready to construct the three-scale factorization theorem by 
combining Eqs. (§) and (^. Start with the nonleptonic heavy meson decay 
amplitude up to 0{as) without integrating out the W boson. The IR sensitive 
functions are hrst factorized according to Fig. 3(a), such that the diagrams in 
the hrst parentheses are characterized by momenta larger than the IR cutoh. 
Employing Fig. 2(b) to separate Hr, we arrive at the factorization formula 
described by Fig. 3(b). The diagrams in the last parentheses are identihed as 
the hard decay subamplitude H. It is obvious that its anomalous dimension 
is given by 'Jh = —l(i> — iHr- Applying the RG analysis to each convolution 







factor, we derive 


n) = c(t)H(t, t)(j){b, 1/6) exp 



where the Wilson coefficient c{t) represents the exponential in Eq. (j^). The 
cancellation of the p dependences among the three convolution factors is ex¬ 
plicit. The two-stage evolutions from 1/6 to f and from t to Mw have been 
established. We emphasize that the Wilson coefficient appears as a con¬ 
volution factor of the three-scale factorization formula, which is, however, a 
constant coefficient (once its argument /i is set to a value) in the conventional 
approach of effective held theory. 

In the leading logarithmic approximation ci _2 are given, in terms of the 
combination c±(/r) = Ci(/i) ± C2(/i), by 


c±(h) 


as{Mw) 

(h) 


33 —2n j! 




( 11 ) 


with the constants 27 _|_ = — 7 _ = —2, and Uf the number of active quark 
havors. Below we shall employ the more complicated two-loop expressions 
of Cl 2 presented in the Appendix A, that include next-to-leading logarithms 

0 -’ 

At last, we explain how to incorporate the Sudakov factor into the above 
factorization formula. The RG solution in Eq, sums only the single log¬ 
arithms contained in the meson wave function (j). In fact, there exist also 
double logarithms coming from the overlap of collinear and soft divergences. 
Hence, an extra large scale P, the meson momentum, should be added into 
0 as an argument. The scale P is associated with the collinear divergences, 
while the small scale 1/6 is associated with the soft divergences as stated 
at the beginning of this section. Before reaching Eq. (^, one performs the 
resummation for these double logarithms, and obtain 


0(P, 6, /i) = 0(6, /i) exp[-s(P, 6)] . (12) 


e~^ is the Sudakov factor, which exhibits a strong suppression in the large 6 re¬ 
gion. The single-scale wave function 0(6, /i) discussed above is then identihed 
as the initial condition of the resummation for the two-scale wave function 
0(P, 6, /i). For the detailed derivation of Eq. (0), refer to [^, . 
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In summary, the large logarithms \n{Mw/t) are grouped into the Wil¬ 
son coefficients ci^ 2 , and ln(t6) are organized by the resummation technique 
and by the RG method. Combining Eqs. (0) and (0), we derive the hnal 
expression of the three-scale factorization formula 


Hr{Mw, /i)0(x, P, b, n) 


c(t)H(t, t)(j){x, b, 1/6) 


X exp 


-s{P,b) 



where the momentum fraction x associated with a valence quark of the meson 
has been inserted. 


III. The B D^*^7r{p) Decays 


In the conventional BSW approach the branching ratios of the exclusive 
nonleptonic heavy meson decays are parametrized only by the factorizable 
contributions from the external and internal W emissions as stated in the 
Introduction. The associated nonfactorizable contributions, which can not be 
expressed in terms of hadronic form factors, are ignored. The nonspectator 
contributions from the IT-exchange (or annihilation) diagrams, which may be 
factorizable or nonfactorizable, are not included either. However, the naive 
PQCD analysis based on the full Hamiltonian in Eq. (|l|) has shown that the 
nonfactorizable contributions are comparable to the factorizable ones for the 
internal W emissions and for the W exchanges p]. 

In this section we shall investigate the importance of the nonfactorizable 
and nonspectator contributions to exclusive nonleptonic heavy meson decays 
employing the more sophiscated three-scale PQCD factorization theorem de¬ 
veloped in the previous section |^. We evaluate the branching ratios of 
the B —> 7r(p) decays, taking into account the factorizable, nonfactoriz¬ 

able and nonspectator contributions, and letting the Wilson coefficients ci ,2 
evolve according to effective held theory. In this framework the external 
W emissions also give nonfactorizable contributions. The relevant effective 
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Hamiltonian is given by 


ftft = ^V,tV;,lcl^^L)0, + c,(p)02] , (14) 

with the four-fermion operators Oi = {du){ch) and O 2 = {cu)[dh). The full 
Hamiltonian H is then a special case with the choice Ci = 1 and C 2 = 0. 

We hrst study the B —> decays. The analysis of the B —>• 

decays is similar. The factorizable external IT-emission amplitudes dehne the 
B —> D^*'^ transition form factors ^ through the hadronic matrix elements, 

{D{P2)\V^\B{P,)) = P v^Y + - v^Y] , 

{D\P2)\A^\B{P,)) = ^MBMB*[UAd){ri + 1 ) 6 *^ - UAvY • 

-Uii.'ny ■ viv^] ■ ( 15 ) 

Pi (P 2 ), Mb {Mbm) and ui (^ 2 ) are the momentum, the mass, and the veloc¬ 
ity of the B meson, satisfying the relation Pi = M^Ui (P 2 = Mb,{*)V 2 )- 

e* is the polarization vector of the D* meson. The velocity transfer ui • V 2 in 
two-body nonleptonic decays takes the maximal value r] = (1 -|-r^)/( 2 r) with 
r = Mb)(*)/Mb- In the rest frame of the B meson Pi and P 2 are expressed as 
Pi = (Mb/a/ 2) (1,1, Or) and P 2 = (Mb/a/ 2) (1, r^. Or) [0]. For the analysis 
below, we dehne ki (/C 2 ) the momentum of the light valence quark in the B 
(pF)) meson. ki may have a minus component fc/, giving the momentum 
fraction xi = /cf/Ph, and small transverse components kir. /c 2 may have a 
large plus component fc/, giving X 2 = k^ jP^ 1 s-nd small k 2 r- The pion then 
carries the momentum P 3 = Pi — P 2 , whose nonvanishing component is only 
P 3 ". One of its valence quark carries the fractional momentum a^aPs, and 
small transverse momenta ksr. In the inhnite mass limit of Mb and Mb(*) 
the form factors with i = -T, —, V, Ai, A 2 , and A^ obey the relations 


'b-l- 'bj4i Cas 


y = u = 0. (16) 

which is normalized to unity 


^ is the so-called Isgur-Wise (IW) function 
at zero recoil 77 —>■ 1 by heavy quark symmetry. 

include the contributions from the hadronic matrix element of Oi 
shown in Fig. 4(a) and from the color-suppressed matrix element of O 2 in 
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Fig. 4(b). Therefore, their factorization formulas involve the Wilson coef- 
hcient Oi = Ci + C 2 /N. We dehne the form factors for the internal 
IF-emission diagrams, which include the factorizable contributions from the 
matrix elements of O2 in Fig. 4(c), and from the color-suppressed matrix 
element of Oi in Fig. 4(d). These form factors then contain the Wilson co¬ 
efficient a2 = C2 + Ci/N. Similarly, we dehne the form factors for the 
VF-exchange diagrams, which include the factorizable contributions from the 
matrix elements of O2 in Fig. 4(e), and from the color-suppressed matrix 
element of Oi in Fig. 4(f). Hence, these form factors also contain the Wilson 
coefficient 02 . 

For the nonfactorizable contributions to the B decays, the pos¬ 

sible diagrams are exhibited in Fig. 5, which correspond to those in Fig. 4. 
Figs. 5(a), 5(c), and 5(e) do not contribute at 0{as) simply because a trace 
of odd number of color matrices vanishes. Hence, all the nonfactorizable 
contributions come from Figs. 5(b), 5(d), and 5(f), denoted by the ampli¬ 
tudes •^d*\ and ■M^f\ respectively, and are thus color-suppressed. 

Their expressions are more complicated, and can not be written in terms of 
hadronic form factors. The amplitudes for the nonfactorizable external 
W emissions depend on the Wilson coefficient C2/N. They have the same 
expressions for the charged and neutral B meson decays, because replacing 
the spectator u quark in the B~ meson by the d quark does not change the 
Feynman rules. The amplitudes for the nonfactorizable internal W 

emissions contain the Wilson coefficient Ci/N. Af for the nonfactorizable 
W exchanges involve the Wilson coefficient Ci/N. 

The decay rates oi B D^*^7r have the expression 

r* = , (17) 

l/OTT r 

where i = 1, 2, 3 and 4 denote the modes B~ D^n~, B^ D^n~, 
B~ —> D*^7r~ and B^ D*~^7r~, respectively. With the above form factors 
and the nonfactorizable amplitudes, the decay amplitudes Afi are written as 

All = — {1 — r)^-] + forint +-Mb +-Md , (18) 

AI 2 = fn[{i + — {I — r)^^] + fs^exc + -Mb + -Mf , (19) 

AI3 = + ^^)^Ai - ( 1 -r)«A2 T^Ag)] 
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( 20 ) 


+/C-S, + Ml + M’a, 

Mi = i^/,|(l + >-)5a, - (1 

+ Me*xc+Xj* + Ai;, (21) 


where /b, and are the B meson, D^*'> meson, and pion decay con¬ 

stants. We have made explicit that the charged B meson decays B~ —> 
£)Wo 

TT contain the external and internal hh-emission contributions, and the 
neutral B meson decays B^ contain the external hh-emission and 

hh-exchange contributions. 

In the considered maximal recoil region with P 2 S> M^(») /\/2 S> Pp, we 
regard the meson as a light meson for simplicity fl^. Double logarithms 


contained in the B meson, meson and pion wave functions cfys, 0£)(*), 
and respectively, are organized into the corresponding Sudakov factors 
using the resummation technique O, |T^, ^ : 


0B(xi,Pi,6i,/i) = 0B(xi,6i,l/6i)exp[-S'B(/i)] , 
0 ^(.)(a: 2 ,P 2 ,& 2 ,/i) = 0^(.) ( 0 : 2 , 62 , 1 /^ 2 ) exp[-S'£,(*) (/x)] , 

07 r(a: 3 ,P 3 , 63 ,/i) = 0 ,r(a: 3 , 63 , l/ 63 )exp[-S'^(p)] , (22) 


with the exponents 

= s{xiP^,bi)+ 2 f ^7(as(/2)) , 

Jl/bi fl 

5'D(*)(h) = s(o: 2 P 2 "^, 62 )+ s((l- o; 2 )P 2 ’*’)^ 2 )+ 2 /" ^ 7 (q;^(/ 2 )) , 

01/62 

SM = s{x3P^, 63) + s((l - X 3 )P 3 , 63) + 2 / — 7(a,(/2)) . (23) 

Jl/bs fl 

The quark anomalous dimension 7 = —as/Tr, is related to 7 ^ = 27 introduced 
in Sec. II. The spatial extents 6* of the mesons are the Fourier conjugate vari¬ 
ables of kiT- We shall neglect the intrinsic b dependence of the wave functions, 
denoted by the argument b, and the 0 (q;s( 1 / 6 )) corrections, denoted by the 
argument 1/6. That is, we assume 0(x, 6 ,1/6) = (j){x). The wave functions 
4>i{x), i = B, D, D*, and tt, satisfy the normalization 



fi 

2Vq ■ 


(24) 
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The exponent s is written as 
s{Q,b)= [ 


ii/b /i 


In 


Q\ 

fi) 




( 25 ) 


where the anomalous dimensions A to two loops and B to one loop are given 
by 


A = 


B = 


^ Cts 

Cf -h 

TT 

-^ln[ 

3 TT \ 


'67 

TT^ 

10 

9 

3 

27 


Hf + ^/3i In 


37b 


'g 27 B-l' 


(26) 


with Cp = 4/3 the color factor and the Euler constant. The two-loop 
expression of the running coupling constant, 


as(/i) 4 16/3i lnln(/i^/A^) 

TT /3oln(pV^^) /^o ln^(/i^/A^) 

will be substituted into Eq. (]25|), with the coefficients 


(27) 


Po — 


33 — 2n 


/ 


Pi = 


153 - 19nf 
6 


(28) 


and the QCD scale A = Aqcd- 

Combined with the evolution of the hard decay subamplitudes H, the 


variables /i in Eq. (^) are replaced by the hard scales t as shown in Eq. 
leading to the RG invariant Sudakov exponents Ssit), and 

Since large logarithms have been organized, we compute H to lowest order 
by sandwiching Figs. 4 and 5 with the matrix structures (/Pi -f Mb)'^^/V2N 
from the initial B meson, with 75 (/p 2 -|- Mp)/y/2N, /e{/P 2 -|- Mp*)/V2N, 
and 75 /p 3 /-\/ 2 iV from the final D meson, D* meson, and pion, respectively. 

The expressions of all the form factors and nonfactorizable amplitudes 
for the B ^ TT decays are listed below. The form factors i = - 1 -, Ai 
and A 3 , and ^j, j = — and A 2 , derived from Figs. 4(a) and (b), are given by 

/* 1 /*C30 

= 167rCp\/rM^ / dxidx2 / bidbib2db2(pBixi)(j)£)(*){x2)ai{t) 

Jo Jo 

XQ!^(f)[(l CiX 2 r)h{xi,X 2 , bi, 62 , m) + {r + (■Xi)h{x 2 , Xi, 62 , ^i)] 
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( 29 ) 


X exp|-SB(() - Sdm(*)I . 

P 1 POO 

= WuCFy/rM"^ / dxidx2 / bidbib2db2(t)B{xi)4>F){*){x2)ai{t) 
Jo Jo 

xas{t)[CjX 2 rh{xi, X2, &i, ^2) + CjXih{x 2 , Xi, 62, &i)] 

X expl-Ssit) - S'^(*)(t)] , 

with the constants 0 


( 30 ) 


c+ = C = 2 


3 

^“2 + 


/ h -1 

■/7 + 1 


"-2 


c- = -Cl = -2 


^(-2 + 


/h + 1 


^-2 


2 ( 77 + 1 ) ’ 

a. = 0, Cl = -1 - 


2(77 + 1 ) 


77 


Cx 43 „ 


77 


20 P 


- 1 
, ci = 




( 31 ) 


Here 77 takes the maximal velocity transfer 77 = (1 + r^)/( 2 r) as stated before. 
For the behavior of the above form factors at other values of velocity transfer, 
refer to . 

The form factors Cint from Figs. 4 (c) and 4 (d), and Clc from Figs. 4 (e) 
and 4 (f) are written as 

Cil = IQttCfVt-MI [ dxidx3 [ 6id6i&3C?&30i?(a;i)07r(ai3)a2(frnt) 

Jo Jo 

xasitint) [(1 + 013(1 -r^))hint(iri, 0:3, 61,63, TTlint) 
+Ci2a;ir^hint(a;3, ici, 63, 61,7771^) exp(tint) - 5 '^(frnt)], ( 32 ) 

p1 poo 

H = 167rCFi/rM| / dX2dX3 / b2db2b3db3(j)F(*){x2)(j)n{x3)a2{texc) 

Jo Jo 

xas(texc) (013(1 - rl - ClH^exc(oi2, OI3, 62, 63, 777 exc) 
+0l2hexc(0l3, 012,63, 62, 777 exc)]exp[-S'ji(*)(fexc) “ 5 ',r( 4 xc)] , ( 33 ) 


with the constants Cint = Cexc = —Cint = “Cxc = 1 - the derivation of Ci 


(*) 

int 


we have assumed that ki has a plus component kf = XiP^. It is obvious that 
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"^int and are exactly the 5 —>■ tt and Dd) —>. 71 - transition form factors, 
respectively, evalnated at maximal recoil. 

In Eqs. (p9D, (p0|), (|3^) and (Q) the fnnctions h’s, obtained from the 
Fourier transform of the lowest-order H, are given by 


h{xi,X2,bi,b2,m) 


hint (0:1,0:3,61, 63, mint) 


hexc {x2,X3,b2,b3,m 

exc) 


Ko {y/xiX2mbi) 

X [6'(6i - b2)Ko {^/x2mbl) Iq {y/x^b2) 

+e{b2 - bi)Ko (y^imfea) lo iV^2mbi )], (34) 
h(o:i,o:3,6i,63,mint) , (35) 


TT 


{y/X2X3mey,cb2) 

X 0(62 - b3)H^^^ {V^3'me^cb2) Jo iv^sme^cbs) 

+0(63 - b2)H^^^ {^yx3me xc 63) Jo {\/X 3 me xcbz 


( 36 ) 


with m = and mint = mexc = (1 ~ We observe that the W- 

exchange contributions are complex due to the exchange of time-like hard 
gluons. The hard scales t are chosen as 


ma.x[^/xlm, ^/x 2 m, 1/bi, I/62) 

( 37 ) 

max{y/x/m~t, ^/x^^^t, I/61,1/63) 

( 38 ) 

max{y/x^m^, a/T^uw, I/62,1/63) • 

( 39 ) 


where we consider only the energies and ^Jxjfn from the longitudinal 

momenta of the internal quarks in the diagrams of Fig. 4, because the gluon 
energies yJXiXjim are always smaller. The scales 1/bi are associated with 
the transverse momenta. Then Sudakov suppression guarantees that the 
main contributions come from the large t region, where the running coupling 
constant as{t) is small, and perturbation theory is reliable. 

For the nonfactorizable amplitudes, the factorization formulas involve 
the kinematic variables of all the three mesons, and the Sudakov exponent 
is given by S' = + Sdm + S'^. The integration over 63 can be performed 

trivially, leading to 63 = 61 or 63 = 62- Their expressions are 

mI*'^ = 327 rA/^CirA/rM| [ [dx] [ 6id6i62d620B(a:i)0£){*)(a:2)07r(a:3) 

Jo Jo 
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X { i^b'’ ) exp [- S ) I b3=fe2 ] 

x[(l - r^)(l - X3) -Xi + ci*\r - r^){xi - X2)]h^b\xi,bi) 

- ^ exp [- ^ (tf 4 163=62 ] 

X [(2 - r)xi - (1 - r)x 2 - (1 - r^)x 3 ]h^^\xi, 6i)| , (40) 

= 327rA/^Cir^/rM| [ [dx] [ &ic?&i& 2 C?& 20 B(a;i) 0 £)(*)(a; 2 ) 07 r(a; 3 ) 

Jo Jo 

X Cd'’ { ^ exp [-^) I b 3 =fe 3 ] 

x[a;i -X2 - X3{1 - r'^)]h^J\xi,bi) 

+«* ^ exp[-^(4^^)Ib 3 =fei] 

x[(a:i +a;2)(l + 4*^4 “ ^]hT{xiM)^ , (41) 

= 3271 V^C F ^/r J^[dx] bidbib2db2(l)B{xi)(j)D(*){^2)4>iT{x3) 

^ exp [-^(t 44 163=62] 

X [xi(l + C/*^r4 - C/*^a:2r’^ - X3(l - r^)]hf\xi, h) 

^ exp[-^(4^4lfe3=fe2] 

X [(xi + X2)(l + C/*V4 - Cf^r‘^]hf{xi, 6i)| , (42) 

with the definition [dx] = dxidx 2 dx 3 . The constants are Cb,dj = ~Cb,d,f ~ 
The functions j = 1 and 2, appearing in Eqs. (^0l)-(^^, are written 
as 

= [e{b3-b2)Ko{BMBb3)Io{BMBb2) 

+ 0(62 - bi)Ko {BMBb 2 ) lo (BMBbi)] 
f Ko{B,Msb 2 ) fori?,> 0 \ 

^-fHl^\\B,\MBb 2 ) forS, <0 y ^ ^ 

hf = [e{b3-b2)Ko{DMBbi)h{DMBb2) 
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+9{b2 - bi)Ko {DMBb2) h {DMBbi)] 
/ Ko(DjMBb 2 ) for B, > 0 \ 

"" [ fi7«(|D,|MB62) forD, <0 J ’ 

hf = z^f9(bi-b2)Ili^U-P'^Bbi)Jo(FMBb2) 

+ 9 (b 2 - (FMBb2) Jo (FMBb,) 

/ Ko{FjMBbi) for Fj > 0 \ 

^ fff«(|F,|MB6i) forF, <oJ ’ 

with the variables 


( 44 ) 


( 45 ) 


B"^ = XiX2{l - r"^) 

Bl = {xi - X2)xo{l - r"^) F XiX2{l + r"^) , 

Bl = a;ia;2(l + r^) - (xi-a;2)(l-a;3)(l-r^) , 

_D2 ^ F"^ = xiXo{l — r"^) ■, 

Dl = Fl = {xi - X2)xo{l - r"^) , 

Dl = (xi + a;2)r^ - (1 - xi - a;2)a;3(l - r^) , 

F2 = {xi F X2) + {I - Xi - X2)xo{l - r'^) . (46) 

The scales are chosen as 

= max(FMB, |Fj|Mb, 1/61,1/62) , 

= max(F)M b, |-Dj|Mb, 1 / 61 , 1 / 62 ) 

= max(FMB, |F,|Mb, 1/61,1/62) . ( 47) 

Here we include also the gluon energies BMb, DMb, and FMb except for 
the energies |-Bj|Mb, |Fj|Mb, and |F) |Mb of the internal quarks, because 
the former may not be smaller than the latter. 

The corresponding form factors and the nonfactorizable amplitudes for 
the B p decays can be computed in a similar way, and their expressions 

are presented in the Appendix B. The only differences are the matrix struc¬ 
tures of the p meson in the calculation of the hard decay subamplitudes, and 
the extra transverse-mode contributions from the pt meson, except for the 
longitudinal-mode contributions from the pi meson. The matrix structures 
associated with the hnal pi and pt mesons are /P^/\/2N and /e /P 3 / \J2N 
with e ■ P 3 = 0, respectively. We also assume a vanishing p meson mass. 
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In the evaluation of the various form factors and amplitudes, we adopt 
Gp = 1.16639 X 10“® GeV“^, the decay constants fp = 200 MeV, fp = 
fp* = 220 MeV, /tt = 132 MeV H, and /„ = 200 MeV, the CKM matrix 


elements \Vcb\ = 0.043 |T^, |T^ and \Vud\ = 0.974, the masses Mp = 5.28 GeV, 
Mp) = 1.87 GeV, and Mp,* = 2.01 GeV [^, and the {B~) meson lifetime 
Tpo = 1.53 (tb- = 1 .68) ps 1^. As to the wave functions, we employ the 
model 

, . X _ Ns p)M _ 3^(1 - x)'^ _ 

^ 167 r 2 + Gb ^ b (*)(1 - a ;) ’ 


(48) 


for the B and mesons, and the Ghernyak-Zhitnitsky models from QGD 
sum rules p^ . 


4>n{x) 

2 ’ 

(49) 


= ^/p3;(l-3;)[0.25(l-2x)2 + 0.15] , 

(50) 


\2 

= 2 • 

(51) 


for the pion and the p meson, respectively. 

The normalization constant Np and the shape parameter Cp are deter¬ 
mined by two constraints from the relativistic constituent quark model [|T^ . 
They are given by Np = 604.332 GeV^ and Cp = —27.5 GeV^, which cor¬ 
respond to fp = 200 MeV listed above. The shape parameters Cp and Cp 
are adjusted such that our predictions for the branching ratios of the various 
modes of i? ^ fall into the errors of the experimental data 0 shown in 

Table I. We determine Cp = —3.372 GeV^ and Cp* = —3.772 GeV^, and the 
corresponding normalization constants Np = 92.85 GeV^ and Np* = 119.51 
GeV^ from fp = fp* = 220 GeV. Results along with those from the naive 
PQGD formalism based on the full Hamiltonian H [ie. with ci = 1 and 
C 2 = 0), are exhibited in Table I. We hnd that the predictions for B meson 
decays from these two approaches are close to each other. For comparision, 
we quote the BSW results from [|, (BSWI), and from 0 (BSWII), in 
which a modihed pole ansatz is employed to make the extraction of the ratio 
02/01 less mode dependent. 

Different kinds of contributions to the decay amplitudes ATi in Eqs. (|T^- 
(pip associated with the decays B —> are presented in Table II. It is 
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f*) 

obvious that the nonfactorizable internal Ph-emission amplitudes play 
an essential role for the explanation of the branching ratios of the B 
Di*) 7T decays: In the charged B meson decays is about 20% of the 

factorizable external Ph-emission contributions, while in the neutral B meson 
decays only the factorizable external PP^ emissions dominate, and all other 
kinds of contributions are small. Hence, the branching ratios of the former 
are predicted to be (1.2)^ x {tb- /tbo) ~ 1.6 times of the latter, which is well 
consistent with the data. This conclusion differs from the previous one drawn 
in §, which is based on the naive PQCD formalism: The factorizable internal 
W emissions give 20% of the factorizable external Ph-emission contributions, 
and are responsible for the ratio of the charged B to neutral B decay rates. 
Therefore, we emphasize that though the values in Column I and in Column 
II of Table I are close, the relative weights of the various contributions change. 
The IT-exchange contributions are always negligible, which are only about 
5% of the external IT-emission amplitudes. If the conventional factorization 
hypothesis for nonleptonic B meson decays is correct, only the diagrams in 
Fig. 4 are considered. However, our analysis has indicated that Figs. 4(c) and 
4(d) give small contributions. This is the reason the naive choice /i = for 
the arguments of the Wilson coefficients in the BSW model can not match 
the data. 

The results for the branching ratios of the B —>■ D^*^p decays listed in 
Table I are also satisfactory. Note that after fixing the meson wave 
function from the data of B ^ there is no free parameter left in the 

analysis of the B D^*^p decays. Hence, the consistency of our predictions 
with the data is nontrivial. The scale dependence of the modified PQCD 
formalism can be tested simply by substituting 2t for t in the factorization 
formulas. The predictions decrease a bit as shown in Table I. In the con¬ 
ventional approach of effective field theory the substitution of Mh by 2Mb 
for the arguments of the Wilson coefficients ci ^2 results in a more than 20% 
difference. Hence, the scale setting ambiguity is indeed moderated in our 
approach. 


IV. The D Decays 
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We have stated that the naive choice of the BSW parameters ai = 
Cl (Me) + C 2 {Mc)/N and 02 = C 2 (Mc) + ci{Mc)/N can not explain the data 
of charm decays. To do it, the large N ansatz ai = Ci(Mc) ~ 1.26 and 
02 = C 2 (Mc) ~ —0.51 must be assumed [jl[. However, the same ansatz 
02 = C 2 {Mh) ~ 0.11 does not work for bottom decays, because the best £t 
to the experimental data gives 02 ~ 0.22 |^. We argue that the large N 
ansatz is the consequence of the factorization hypothesis employed in the 
BSW model. If the nonfactorizable contributions along with the evolution of 
the Wilson coefficients are taken into account, such an ansatz is not necessary. 
In this section we apply the three-scale factorization theorem to the decays 
D TT, and explore in details how the contributions from Figs. 4 and 5 

vary, when they are evaluated at different energy scales. In our approach the 
factorizable and nonfactorizable contributions change with the characteristic 
scales t of the decay processes in different ways, such that their combination 
can explain both the bottom and charm data. That is, our work provides a 
unihed viewpoint to the exclusive nonleptonic heavy meson decays. 

Before proceeding with the calculation of the decay rates, we emphasize 
that the applicability of PQCD to D meson decays with Md = 1.87 GeV is 
marginal. It has been shown that the PQCD analysis of exclusive processes 
with the Sudakov effects included is reliable for the energy scale above 2 GeV 
2% . Therefore, we concentrate only on the mechanism of the destructive 


interference involved in charm decays. For this purpose, it is enough to 
consider the ratios of the charged D meson decay rates to the neutral D 
meson decay rates. 


Ri = 


B{D- K^tt-) 
K+n-) 


R2 = 


B{D- K*^7i-) 
B{D^ K*+7i-) ’ 


(62) 


instead of the individual branching ratios. 

The D TT decays occur through a similar effective Hamiltonian 




(53) 


where the four-fermion operators are Oi = {du){sc) and O 2 = {su){dc). 
The analysis of the nonleptonic B meson decays in the previous section can 
be copied to the D meson decays directly. The expressions of the decay 
rates Fj are similar to Eq. (Q. but with the subscripts i = 1, 2, 3, and 4 
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denoting the modes D~ —^ K^Tr~, > K^tt~, D~ and Z)° —>• 

, respectively. At the same time, the CKM matrix element |Ksl = 1.01 
is substituted for | 14 b|, and the meson masses and for Mb and 

Mb(»), respectively. In all the form factors and nonfactorizable amplitudes 
the kinematic variables of the B meson are replaced by those of the D 

meson. The Z)° {D~) meson lifetime is Tbo = 0.415 {td- = 1.05) ps 
[^. The D meson wave function has been determined in the study of the 
B meson decays. For the kaon, we have the masses Mk = 0.497 GeV and 
Mk* = 0.892 GeV, the decay constants fK = 160 and fx* = 220 MeV, and 
the model wave functions derived from QGD sum rules p^. 



= ^/^x(l-x)[3.0(1-2 x)2 +0.4] , 

(54) 

<pK*ix) 

= ^/ic*a;(l - x) [0.5(1 - 2x)^ + 0.9] , 

(55) 

(/>k4x) 

= \/6fx*x(l — x)[0.7 — (1 — 2x)^] , 

(56) 


for the K, and mesons, respectively. Note that we take cj)^* as the 
K* meson wave functions throughout the analysis of the D meson decays for 
simplicity. Then all the factorization formulas in Sec. Ill can be adopted 
directly without further modihcation. Gompared to the pion wave funtion 
(I)t, in Eq. (^91), (pxs do not possess dips at the middle of the momentum 
fraction x. 0^* has a single hump at x = 1/2, differing from the behavior of 
(j)K and 

Because of the smaller D meson mass, the transverse degrees of freedom 
are more important in the dehnitions of the hard scales t. Hence, we choose 
the maximum of the scales 1/bi for the arguments t of the Wilson coefficients. 
In this case Sudakov suppression is weaker, and thus insufficient to diminish 
the contributions from the region with t close to Aqcd, where ci ,2 diverge. 
To have meaningful predictions, a lower bound tc = {1 + e)AQCD must be 
introdiced for the variables t in the numerical analysis, where e is a small 
number. We then have t = max(l/6j,G)- The results of Ri and R 2 for 
e ~ 0.0002 are exhibited in Table III, which are well consistent with the 
data. Note that tc can be regarded as one and the only one free parameter 
in the analysis of the D meson decays. Therefore, the simultaneous £t to Ri 
and i ?2 is not trivial. 

The contributions from different diagrams are listed in Table IV. The W- 
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exchange contributions to the neutral D meson decays are negligible as in 
the neutral B meson decays. It is easy to observe that with t running to be¬ 
low Me, the factorizable internal-hh emission contributions to the charged D 
meson decays become very negative due to the evolution of 02 , and overcome 
the positive nonfactorizable internal hh-emission amplitudes Note that 

the factorizable internal-iy emission contributions are positive and small in 
the B meson decays. The naive PQCD formalism based on the full Hamil¬ 
tonian H without considering the evolution of the Wilson coefficients |^, § 
can not account for this sign change, since the corresponding coefficient 02 
is always equal to 1/N. It then predicts that the charged D meson decay 
rates are larger than the neutral ones (the values of Ri and R 2 in Column I 
of Table III are greater than unity) as in the B meson case. 


V. The B Decays 


As mentioned in the Introduction, it has been very difficult to explain 
the ratios R and Rl associated with the B —>• decays simultane¬ 

ously, which were dehned in Eq. (^), in the BSW framework based on the 
factorization hypothesis. We have observed in Secs. Ill and IV that the non¬ 
factorizable contributions play an essential role in the decays B D^*'>Ti[p) 
and D 71. Therefore, it is expected that the nonfactorizable effects 

are also important in the decays B —>■ In fact, it has been sus¬ 

pected that the discrepancy between model-dependent BSW predictions and 
the experimental data 0 is attributed to the breakdown of the factorization 
hypothesis . 

In this section we apply the three-scale factorization theorem to the B 
decays, and show that our predictions for the branching ratios of the 
various decay modes and for R and Rl are in good agreement with the data. 
Similarly, the decays B —> J/ip occur through the effective Hamiltonian, 


= ^v'dV;;ici(fi)o, + c2(,^)02], 


(57) 


with the four-fermion operators Oi = {sc){cb) and O 2 = {cc){sb). The rele- 
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vant decay rates have the expression 

r. = , (ss) 

with r = Mj/^/Mb, Mj/^ = 3.096 GeV being the J/'ifj meson mass. The 
subscript i = 1 denotes the modes B~ J/'tpK~ and —>■ J/'ipK^, which 

possess the same factorization formulas, and i = 2 denotes B~ J/iIjK*~ 
and B'^ —>■ J/ipK*^. Since the decay amplitudes M.i contain only the internal 
hh-emission contributions from Figs. 4(c) and 4(d) as the factorizable part, 
and from Fig. 5(d) as the nonfactorizable part, their expressions are given 
by 


Ml 


(59) 

-M 2 


(60) 

Ml 

= fMM + , 

(61) 


where the superscripts L and T denote the logitudinal and transverse modes, 
B —> and B —> respectively, and fj/^ = 390 MeV is the J/il) 

meson decay constant 0. Eq. (|59D is similar to Eq. (^) and Eqs. (^) and 
(1^ to Eq. (^Of) , but with the external hF-emission contributions dropped. 
The type of the mode B^ —»• J/ipK^ corresponds to that of B^ —> 
which was not considered in Sec. III. Note that the B^ —> D^ir^ decay 
involves not only the internal PF-emission but the PF-exchange diagrams. 

Employing the matrix structure /^(/p 2 + Mj/^)/^/2N for the hnal J/ijj 
meson, and the vanishing kaon masses = Mk» = 0 for simplicity, the fac¬ 
torization formulas of the form factors and of the nonfactorizable amplitudes 
are derived straightforwardly. They are written as 


>int 


_AJ/i>)L 

Sint 

lQ^^CF^/rMB 


dxidxs 


AJIi^)T 

Sint 


/ bidbibsdbs(j)B{xi)(j)^(.){x3)a2{tint 

JV Jo 

XCrs(tint ) 1^(1 -P 0:3(1 - r‘^))hint{xi, X3, bi, 63, mint) 

-Xir‘^hintix3,Xi,b3,bi,mint)] exp[-5B(fint) - ^^(tint)], 

p1 poo 

= 327rCi7’V^M| / dxidx3 / 5id5i53C?&30s(a;i)0^*(a;3)a2(hnt) 


(62) 


xasitint)rhint{xi, X3,bi,b3,mint) exp[-SBitint) - 5'i^(tint)], (63) 
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= -xf . 

X {a*(tf ^ ^ exp[-^( 4 ^^) Ib 3 =bi] 
x[2 -r^ -2(a;i + X 2)(1 

/'j-(2)', 

+«s exp [- S {tf) I fe3=b, ] 

x[4xi - - 2 x 2(1 + r4 - 2 x 3(1 - r‘^)]hf'{xi, 6i)| , (64) 

= 327 iV^CFVrMl [dx] 6i(i6i62(^&20B(a;i)0j/p(a;2)0K*(2:3) 

X{as(trf^) ^ exp [- 5 )Ib3=b^] 
x2r(l - xi - X2)/ii^43^*>^*) 

-as (tf ^ ^ exp [-^( 4 ^^) 163=61 ] 

x 2 r(l-xi +X2)/6d^4^i)^*)| ) ( 65 ) 

The total Sudakov exponent for the nonfactorizable amplitudes is given by 


S = Sb + Sj/^ + Sk, where Sj/^ (Sk) has the same expression as Sf)(*) {S-^) 
in Eq. (^). The hard functions hint and j = 1 and 2, are the same as 

those appearing Eqs. (^5]) and (^^, but with the arguments 

mint = 

= 

Ml - , 

XiX3(l - r4 , 

(66) 

(67) 

Dl = 

(1 - X 2 )xi(l + r^) - (3 - 2X2 - ^ 2 )^ 



+ (xi +X 2 - l)a:3(l -^‘^) , 

(68) 

Dl = 

XiX 2 (l + r^) + (xi - X 2 )X 3(1 - r^) + (x 2 - -)r^ . 

(69) 

The hard scales t 

are also similar to those in the analysis of the B ■ 

DMjt 


decays but with the insertion of the above arguments. 
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To evaluate the form factors and the nonfactorizable amplitudes, we need 
the information of the J/V’ meson wave function. Unfortunately, there are not 
yet convincing models for them. However, it should be most possible that the 
two charm quarks in the J/ip meson carry equal fractional momenta. Hence, 
we assume, for convenience, that the wave function for the ( meson 

with transverse polarization possesses the same form as oc x^(l — for 
the pt meson in Eq. (0), because has a maximum at x = 1/2. We further 
assume that for the {J/'iP)l meson with longitudinal polarization is 

proportional to x"'(l — x)” with n a free parameter, which will be determined 
by the data of the decay B —>• J/%IjK. That is, we propose the model wave 
functions. 




( 70 ) 

( 71 ) 


The constant Njj^ is related to the normalization condition 

^ dx(t)^j/Yx) = ^ . (72) 

We stress that the particular form of the J/'ip meson wave functions are not 
important. We have tried other models, such as x(l — x)exp[—(1 — 2x)^], 
and found that it works equally well. The kaon wave functions have been 
shown in Eqs. (|5lD-(p6D. 

The experimental data of the branching ratios of the B —> J/ipK^*'^ de¬ 
cays, and of R and Rl |0] are exhibited in Table V. We determine the 
parameter n = 1.25 from the branching ratio B{B —>■ J/'ipK), and then the 
normalization = 0.858 GeV^ from Eq. (^). After hxing the J/'ip meson 
wave functions, we evaluate the branching ratios of the decays B J/'ll) 
and B J/rjjK^. Results and the corresponding factorizable and nonfac¬ 
torizable contributions are presented in Table V and VI, respectively. Ob¬ 
viously, both the branching ratios of the various decay modes, and R and 
Rl are explained successfully. If the nonfactorizable amplitudes are 

excluded, we shall have R = 1.48 and Rl = 0.92, which is too large. It 
implies that the nonfactorizable contributions are indeed essential for the 
decays B —»• 
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V. Discussion 


In this paper we have developed a modihed PQCD formalism for the study 
of the exclusive nonleptonic heavy meson decays, which embodies effective 
field theory and factorization theorems. It involves three scales: the W boson 
mass Mw, the characteristic energy t of the decay processes, and the trans¬ 
verse extent b of the mesons. The evolution of the Wilson coefficients from 
Mw to t and of the Sudakov factor from t to 1/b are established to make the 
factorization formulas explicitly /i-independent. The factorizable, nonfactor- 
izable and nonspectator contributions from the external W emissions, the 
internal W emissions, and the W exchanges are all taken into account, and 
have been evaluated reliably. We emphasize again that the Wilson coefficient 
appears as a convolution factor of the factorization formulas in our analysis, 
instead of a constant coefficient as in the conventional approach of effective 
held theory. 

Basically, the two main controversies in the exclusive nonleptonic heavy 
meson decays, ie. the extraction of ai ,2 from bottom and charm decays, and 
the simultaneous explanation of R and Rl, have been resolved by our for¬ 
malism: The nonfactorizable external hT-emission contributions alone, 
which are 20% of the factorizable one, account for the data of the B —> 
Di*) 7r(p) decays. The factorizable internal W-emission contributions, be¬ 
coming negative enough to overcome in the D meson case, successfully 
explain the data of the D K^*'^7r decays. That is, the evolution of the 
Wilson coefficients can lead to the necessary constructive and destructive 
interferences involved in bottom and charm decays. While it is the nonfac¬ 
torizable contributions that make trivial to account for the ratios R and R^ 
associated with the decays B J/ipK‘^*\ 

Note that the free parameters contained in our formalism are less than 
the decay modes considered. Hence, the match of the theoretical predictions 
with the experimental data is nontrivial, and indicates that we may have 
explored the correct mechanism responsible for the nonleptonic heavy meson 
decays. It is worthwhile to compute other decay modes, whose consistency 
with the data will further justify our approach. Inclusive nonleptonic heavy 
meson decays are another important subject to which our formalism can be 
applied. The nonfactorizable soft corrections U will give the hne tuning of 
our predictions. These topics will be investigated in separate works. 
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A Two-loop expressions of 

In this appendix we present the expressions of the Wilson coefficients ci, 2 (/^) 
to two loops, which are given in terms of c± = ci ± C 2 by ^0 


c±(h) = 


Q^s(h') j 


(/^) 


d± 


1 + _ 4 ) 


(Al) 


with the constants 




7i‘> 


w 

A =F 1 
= ± 12 - ^ 


<7+ = 


(0) 

7 ± 

w 


2N 


( 1 ) NtI 

7 + = - 

2N 

NtI 


57 19 4 

-21 ± — ^ —N ±-nf- 2/?ofc± 
A 3 3 ^ ^ 


= 


2A 


[±11 ± k ±\ 


(A2) 


The scheme dependent parameters k± are 


= 0 NDR 

= t4. HV (A3) 

The constants Po and Pi have been defined in Eq. (p8|). In this paper we adopt 
the NDR scheme. However, we have tested the sensitivity of onr predictions 
to these two schemes, and fonnd that the scheme dependence can be absorbed 
into the meson wave fnnctions. Namely, the wave fnnctions vary with the 
scheme snch that the predictions almost remain the same. 




















When the scale fi in as{fi) evolves from above Mb to below Mb, the 
flavor number Uf changes from 5 to 4. A similar change from nj = 4 to 
3 occurs as /x evolves from above Me to below M^. To make continuous 
at these thresholds, Aqcd must change accordingly. However, again, the 
dependence on Aqcd can also be absorbed into the wave functions, such that 
our predictions are insensitive to whether the continuity conditions of are 
implemented or not. Hence, for simplicity and within the errors of the data, 
we assign the value Aqcd = 0.2 GeV, and n/ = 4 for bottom decays and 
Uf = 3 for charm decays in the numerical analysis. 


B Factorization of the B D^*^p Decays 


The factorization formulas for the B ^ D^*^p decays can be derived straight¬ 
forwardly. The only differences from the B —>■ case are the matrix 

structures of the p meson in the calculation of the hard decay subamplitudes 
H, and the extra contributions from the transverse modes involving the pt 
meson, as stated in Sec. HI. 

The decay rates are given by Eq. (|T^) with i = 1, 2, 3 and 4 representing 
the modes B~ D^p~ , B^ B~ D*^p~ and B^ D*^p~ , 

respectively. The decay amplitudes M.i are written as 


Ml 

~ /p[(l + ~ (1 ~ + fD^int — Mb + Md , 

(Bl) 

M2 

— /p[(l + ^)'^+ ~ (1 ~ + /s'Cexc — Mb + Mf , 

(B2) 

-Ms 

- 2^ /p[(l + r)^Ai (1 +CA3)] 



+fD*^L + Ml-M:, 

(B3) 

CO 

= 

(B4) 

M^ 




+/bCxc + Ml - M} , 

(B5) 

Ml 

= mUc + MJ*, 

(B6) 


where the superscript L (T) denotes the longitudinal (transverse) mode B —> 
D* pBT) _ yy-g have used the the same notations as those for the B —>■ 
decays without confusion. 
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The form factors i = +, —, Ai, A 2 and ^ 3 , related only to the B —> 
D^*'> transitions, are the same as those associated with the B decays. 

The form factors and and the nonfactorizable amplitudes Ai^Xf 
similar to those in the B —> decays, but with the pion wave function 

07 r(a^ 3 ) replaced by the pi meson wave function (j)p{x 3 ) given in Eq. (|50|) . 
The Sudakov exponents Sp for the p meson and S-,^ for the pion are the same. 
Below we give only the form factors and and the nonfactorizable 
amplitudes Af involved in the transverse modes B D*pt, 

p 1 poo 

= 167rCFA/rM| dxidxs bidbib5dbs(j)B{xi)(j)^{x3)a2{tint) 

xasitint)rhint{xi,X 3 , bi, 63, mint) exp[-S's(tint) - S'p(tint)] , (B7) 

pi poo 

= IbTrCiJ’v^Ml dx2dx3 &2d62&3d630D*(a;2)0j(a:3)a2(4xc) 

xas(tey,c)r'^hey,c{x 2 , X3,b2,b3, rrie^c) exp[-SD*(t exc) *S^p(^exc)] •> 

(B 8 ) 

M^* = 3271 [ [dx] [ bidbib 2 db 2 (j)B{xi)(j)D*{x 2 )(j)^{x 3 ) 

JO Jo 

^ exp[-^(ty^)|fc3=fejr[l -xi - X2]h^Xi^ubi) 

+(^s{td^) exp[-S{t^a'^)\b,=b,]r{xi - X2)h‘'a\xi,bi)j , 

(B9) 

M'f* = 3277\/^CF\frM% J [dx] J bidbib2db2(f)B{xi)(j)D*{x2)4X{x3) 

x|as( 4 ^^)—];;^exp[-S'(t^^^)| 63 = 62 ]r^[l -xi - X2]hX\xi,bi) 

^ exp[-5'(t^^^)|b3=b3]r^(a;i - X2)hf {xi,bi)^ , 

(BIO) 

with the Sudakov exponent S = Sb + So* + Sp. The functions hint, and 
\ i = 1 and 2, and the hard scales t have been defined in Sec. III. 
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Table I. Predictions for the branching ratios of the B D^*^Ti{p) decays 
from the PQCD formalism based on H (I), on Hes (II), on Hes bnt with 
the hard scale t replaced by 2t (III), and from the BSW model with the 
parameters Oi = 1.15 and 02 = 0.26 (BSWI) 

02/01 = 0.224 (BSWII) 0. The CLEO data 0 are also shown. 


17 and with oi = 1.012 and 


modes 
B- D^tt- 
130 ^ D+ 7 T- 
B- D*^7r- 

B- D^p- 
5° ^ D+p- 
B- D*^p- 
BO ^ D*+p- 


I(%) 

II(%) 

III(%) 

BSWI(%) 

BSWII(%) 

data(%) 

0.52 

0.50 

0.47 

0.57 

0.51 

0.534 ±0.025 

0.33 

0.33 

0.31 

0.35 

0.28 

0.308 ±0.026 

0.49 

0.48 

0.46 

0.56 

0.56 

0.497 ± 0.044 

0.32 

0.32 

0.30 

0.34 

0.27 

0.304 ±0.024 

1.34 

1.21 

1.16 

1.07 

1.11 

1.022 ±0.067 

0.63 

0.68 

0.62 

0.82 

0.69 

0.861 ±0.078 

1.34 

1.62 

1.33 

1.27 

1.48 

1.444 ±0.134 

0.58 

0.83 

0.69 

0.93 

0.83 

0.844 ±0.071 
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Table II. Contributions to the B decays from the factorizable 

external W emissions and internal W emissions (or W exchanges), and from 
the nonfactorizable amplitudes /■ in Eqs. ([I8|) - (^ID . The unit is 10“^ 
GeV. 


amplitudes 

external W 
(factorizable) 

internal W 
(factorizable) 

<’ 


Ml 

95.1 

2.5 

-5.3 + 14.8* 

18.5 - 10.4i 

Ms 

86.5 

2.6 

6.6 — 20.6* 

17.0-10.8* 

amplitudes 

external W 
(factorizable) 

W exchange 
(factorizable) 


<’ 

M 2 

95.1 

—0.6 + OAi 

-5.3 + 14.8* 

2.2-3.1* 

M 4 

86.5 

-0.6+ 1.3* 

6.6 — 20.6* 

3.0-3.1* 
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Table III. Predictions for the ratios Ri and R 2 associated with 
Ki*) TT decays from the PQCD formalism based on H (I) and on 
The experimental data are also shown. 


the D 
(II). 


modes 

I 

II 

data 

Rl 

4.96 

0.74 

0.72 

R2 

5.00 

0.35 

0.38 
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Table IV. Contributions to the D —> 71 decays from the factorizable 

external W emissions and internal W emissions (or W exchanges), and from 
the nonfactorizable amplitudes The unit is 10“^ GeV. 


amplitudes 

external W 
(factorizable) 

internal W 
(factorizable) 

Mi” 


Ml 

368.0 

-192.4 

-19.7 + 17.1f 

18.5 - 24.4f 

Ms 

752.0 

-576.7 

44.9 - 10.3f 

119.2-39.8f 

amplitudes 

external W 
(factorizable) 

W exchange 
(factorizable) 



M 2 

368.0 

3.6-2.If 

-19.7 + 17.1f 

-10.7 + 24.7f 

M 4 

752.0 

-2.5 -37.8f 

44.9 - 10.3f 

51.6-5.6f 
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Table V. Predictions for the branching ratios of the B —> J decays 
from the PQCD formalism based on H (I) and on iPefr (II), and from the 
BSW model with the parameters Oi = 1.012 and 02/01 = 0.224 (BSWII) [^]. 

0 


The CLEO and CDF data 


are also shown. 



modes 

I(%) 

II(%) 

BSWII(%) 

data(%) 

B- 

^ J/i>K- 

0.11 

0.11 


0.102 ±0.014 

50 


0.10 

0.10 


0.115 ±0.023 

B- 


0.14 

0.15 


0.158 ±0.047 

50 


0.13 

0.14 


0.136 ±0.027 


R 

1.32 

1.47 

1.84 

1.36 ±0.17 ±0.04 

1.32 ± 0.23 ±0.16(CDF) 


Rl 

0.62 

0.56 

0.56 

0.52 ±0.07 ±0.04 

0.65 ± 0.10 ±0.04(CDF) 
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Table VI. Contributions to the B J decays from the factorizable 
internal W emissions and from the nonfactorizable amplitudes The 

unit is 10“^ GeV. 


amplitudes 

internal W 
(factorizable) 


Ml 

126.9 

-30.1 + 4.9f 

Mi 

121.0 

-32.8 + 0.2f 

Ml 

36.7 

34.8 + 30.9f 
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Figure Captions 


Fig. 1. (a) Separation of the infrared and hard 0{as) contributions in 
PQCD. (b) 0{as) factorization into a soft function and a hard decay subam- 
plitude. 

Fig. 2. (a) Separation of the hard and harder 0{as) contributions in effec¬ 
tive held theory, (b) 0 {as) factorization into a ’’harder” function and a hard 
decay subamplitude. 

Fig. 3. (a) 0{as) factorization of a soft function from a full decay amplitude, 
(b) 0 {as) three-scale factorization formula for a decay amplitude. 

Fig. 4. Factorizable external W emissions from (a) the operator Oi and 
from (b) O 2 , factorizable internal W emissions from (c) Oi and from (d) O 2 , 
and factorizable W exchanges from (e) Oi and from (f) 02 - 

Fig. 5. Nonfactorizable external W emissions from (a) the operator Oi and 
from (b) O 2 , nonfactorizable internal W emissions from (c) Oi and from (d) 
O 2 , and nonfactorizable W exchanges from (e) Oi and from (f) 02- 
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Fig.3 

























